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ABSTRACT 



The Center for Naval Analyses Computer War Game, SEALIFT, is a Monte- 
Carlo simulation designed to help study sealift capabilities in an ASW 
environment. A mathematical model of the SEALIFT game is posited to ob- 
tain expected value results approximating those of the SEALIFT game. The 
model is cast in the Fortran terminology of SEALIFT and an effort is made 
to accurately reflect the SEALIFT flow chart logic in its development. 
Comparisons with SEALIFT results to determine the model's reliability and 
accuracy are not made. 
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1. Introduction 



Objective ^ In 1963-1964 the Center for Naval Analyses, CNA, of the 
Franklin Institute constructed a rather large and complex computer war game, 
SEALIFT I ^1 j • This game was used as an analytical tool for the Anti Sub- 
marine Warfare Study, CYCLOPS II. Due to the large computer running time 
of SEALIFT an after the fact need for a "quick and dirty" approximation to 
actual SEALIFT results was felt, and during the summer of 1964 this author 
was privileged to work with the CNA CYCLOPS study group on that approxima- 
tion. This thesis is an outgrowth of that summer* s work. The objective of 
this thesis is to provide the logic for an expected value analysis of that 
portion of the SEALIFT game directed toward calculating transport ship 
losses through submarine action. This thesis is then a "first cut" model 
to approximate a model, and the analysis is executed under simplifying 
assumptions. 

Background. Functionally, the SEALIFT game pits the submarine forces 
of one opponent against the sea transport endeavors of another opponent in 
support of a time and scope limited war overseas. The major game output 
is the number of transport ships sunk during the war and, through inference, 
the effect of these losses on the sealift supply rates. 

Structurally, SEALIFT is an event store Monte-Carlo computer simulation 
and is primarily probabilistic with essentially no kinematics or geometry 
wired in. Although certain tactical procedures were assumed during con- 
struction, they may be circumvented or modified by the input parameter set# 
This set is on the order of four hundred and contains numerous decision 
switches as well as numerical data. It is in the scope of this parameter 
set that one finds the flexibility of the game and, unhappily, the bane of 
forthright analysis. 
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The focal point of SEALIFT is the submarine-convoy encounter and of 
the thirty-two subroutines which comprise SEALIFT, five form the substance 
of this event o The remaining routines are directed toward making the game 
more realistic by introducing other strategies and toward controlling the 
flow of the gameo Figure 1 is a graphical representation of the game® A 
general and more detailed description of the game’s structure and logic is 
quoted from SEALIFT I |^lj in Appendix I« 

Philosophv o All submarine loss events in SEALIFT except those incurred 
in the submarine-convoy encounter are analyzed here utilizing independent 
Bernoulli trial techniques to yield simple expected value resultso To 
effect this result, several event-time dependent routines were omitted from 
consideration and a few were modified to simplify the analysis® 

Preliminary investigation indicated that the submarine-convoy inter- 
action was the most sensitive in the determination of transport ship losses 
and that simple techniques were not sufficiently accurate to determine these 
losses® This encounter has been recast in the framework of a Markov chain 
where each of the sequences of events which may befall a submarine upon 
convoy passage is considered a state of the Markov chain® Primary effort 
has been directed toward accurately effecting the Markov structure, and 
minimal effort toward the accurate reflection of each state within the chain. 
Each state, i®e., sequence of events, is described via independent Bernoulli 
trial terminology, again with approximations® 

This analysis has been made with some trepidation® It may be argued 
that detailed Markov analysis is pretentious in view of the coarse Bernoulli 
approximations; nevertheless, if the results of this expected value analysis 
approximate SEALIFT results, it is hoped this method of analysis will pro«»* 
vide new insight for a more detailed expected value analysis of the entire 
SEALIFT game® 
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FIGURE 1 

GRAPHICAL representation OF THE SEALIFT GAME 













2, Mathematical Preliminaries 



Notational usage in this thesis o The following symbols are used in 
this thesis: 

Symbol Definition 



Pr(A) 


The probability of the event A. 


Pr(A/B) 


The probability of the event A conditioned 
on the occurrence of the event Bo 


E(N) 

H 


The expected value of the random variable No 

The greatest integer less than or equal to a 
The brackets are used, additionally, to en- 
close a matrix when its elements are listed 
in fullo 


Iqy) 


The matrix whose ij th element is q. .* 

ij 


A = B 


A is defined as B, or A is identically equal 
to Bo 


A i B 


A is approximately equal to Bo 


n 


The product, A, x A^ x o*o x A « 
^ 12 n 


E *i 

i=l 


The sum, A. + A_ + o o • + A o 
12 n 


iff 


if and only if 
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Markov chains ♦ It is generally understood that if a physical system 
can exist in a known set of stable states, say Ej^, •••, (finite or 

countable), and if with each state there is associated a probability, 

Pr(E.) = p., that the system is in that state (or goes to that state), then 



the probability that the system transits the states E ., E E . is 

n ^0 -^n 



the product Tf p . • The basic assumption here is transitional indepen- 
k=Q 

dence, i.e,, the chance for transition from a given state to another does 
not depend upon the given state or upon past states* 

Some physical systems cannot be so explained but require the dependence 
of transitional probabilities upon the history of the system* The simplest 
of such dependencies is that the next state in sequence depends upon the 
present or given state only* In this case the outcome, E^, is no longer 
associated with a fixed probability p. , but every pair (E., E.) is associ- 

K 1 j 

ated with a conditional probability, P^jf which denotes the chance that 
the system goes to state E^ given that it was in state E^ at the previous 
step* In general, the transitional probability p_ is a function of the 
sub-subscript in the transit E . , • * * , E . , and the process is called homo- 



geneous if p^^ is independent of this sub-subscript* (Only the homogeneous 

Case will be considered in the expected value analysis and in what follows*) 

If, in addition, the probabilities, p^, of the initial trial i are 

known, the probability that the system transits the states E^ , E^ , ***, E. 

«3 , 



Jq Jjl 



n-1 



is the product p . x Tf P • • • The conditional probability of the 

•^0 k=0 ■^k'^k+1 

n-1 

transit, given the initial state E . , is then Tf P - • 

■^0 k=Q ’^k’^k+1 

The transitional probabilities can be arranged in a matrix, called the 
transition probability matrix, P 
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'1 



(2.D) 



P = 



^11 ^12 ^13 



21 ^22 ^23 •' 



31 ^-^32 ^^33 



where the standard matrix subscripts ij also stand for transition from state 

to state E^.o Together with the initial distribution, p^, the matrix 

completely defines a Markov chain with states E^ , E^ i « « o o It should be 

noted that the elements of each row must sum to one, and that the element 

p. . denotes the chance for transition from state E. to state E. in one step* 
ij 1 3 ^ 

The natural extension is to investigate the concept of transition 



from state E^ to state E^. in more than one step, say in n steps® This 

transition can occur via many paths, and hence the conditional probability, 

p!?\ of finding the system in state E. at the n+r th step, given that it 
^ J 3 

was in state E^ at the r th step, is the sum of the conditional probabilities 
for each possible path® In particular, 

p!V = p.. 



The equation 
( 261 ) 



iJ 

( 2 ) 

ij 

1 

(n) 

Pij 



iJ 



r 



V’ (n-1) 
* ^PirPrj 



and the Chapman-Kolmogorov equation 

I 

P; 



(2o2) (m+n) ^ (m) (n) 

= Lp.. P, 






ix rj 



follow by mathematical induction. 



1 ,,. 



"William Feller, ^ Introduction To Probability Theory And Its Appli- 
cations (New York; John Wiley and Sons, 1957) p. 348. 
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I; 
















The probabilities, can be arranged in a matrix of identical structure 

as that given by equation (2,0) • Coupled with equation (2.1), this suggests 
matrix multiplication - the ij th element of the n step matrix is the inner 
product of the i th row of the one-step matrix with the j th column of the 
(n-l) st step matrix, i.e., 

= P X p'’"^ 



Similarly, equation (2.2) implies 
pm+n ^ pm ^ pn 

Also, if p^ is the probability of the initial state E^, the unconditional 
probability of finding the system in state Ej after n steps is 



(2.3) 



,(n) V 
j “ ^ Pi 



(n) 

X p. . 



A particular n step transition matrix . The expected value analysis 
requires the n step matrix, formed from the (one-step) transition prob- 

ability matrix, P, 



(2.4) 



Pi Pz P3 P4 P5 P6 P7 Pb 



Pi Pz P3 P4 P5 P6 P7 Pb 



00100000 

00010000 

00001000 

00000100 

00000010 

00000001 



= (p..) 



where = p^^ = of the previous notation, etc. If A^, A^, A^ 
are defined such that 
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# 









II 

< 


“Pl P2“ 


’ A^ - 


~P3 P4 


, A3 = 


"P5 P6 


> 

II 


~Pt Pb 




Pi P2_ 




_P3 P4_ 




P5 P6_ 




P7 P0_ 



then 



and 

( 2 . 5 ) 



and in general 



Ai A^ A^ 
0 10 0 
0 0 10 
0 0 0 1 



0 

0 

0 



0 

0 

0 



(Aj^+I ) A2 

I 

0 

0 



n-1 

iY, aJ)A2 

i=0 

I 

0 

0 



(Aj^+DA^ 

0 

I 

0 



(A^+I)A^ 

0 

0 

I 



n-1 

{T a"-)a. 

i=D ^ 



0 

I 

0 



By recursive multiplication, it follows that 



a" = (Pj^ + Aj^, n 5 1 



and therefore 



n-1 . 

(E \ 



i=^0 



0 

0 

I 



n-1 n-1 . n-1 

E ‘ E ‘ ' E <Pi •- P2> 

1=0 1=1 1=1 



i-1 



X + I 



1 - (Pl + P2^ 

X A- + I 

1 - (Pj^ + p^) 
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a 






I 

i 



Letting 



f2(n) = 



1 - (p^ + P 2 ) 

1 “ (Pl + 



n-1 



it follows that 



n-1 . 

I \ 

i=0 



^1 '2 
^1 ^2 



and substitution in (2.4) obtains 



( 2 . 6 ) 



(Pj^ -t- P2) X 





1 - (p, + 


; (n) = 




1 - (pj^ + 


(n) + 1 


P 2 ^ 2 ^"^ 


(n) 


P 2 ^ 2 ^^^ * ^ 


Pi P 2 


f (n) X 


Pi P 2 


1 


^6x2 





P3 P4 P5 P6 P7 P 
P3 P4 P5 P6 P7 P 



'6x6 
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3o The Model for the Expected Value Analysis 

In what follows, a perusal of Appendix I and of SEALIFT I ^ij is 
assumed® 

Methodology ® The heart of any computer game is the flow chart logic, 
and it was via these charts that SEALIFT was submitted to analysis® Prior 
to detailed analysis, the game was examined as a unit to cull those routines 
pertinent to the submarine from those of a bookkeeping nature and other 
routines of non-interest, including those which determined convoy losses 
to non-submarine causes® Of the thirty-two routines which comprise <the 
SEALIFT game, twenty were discarded as irrelevant to the objective of this 
thesis. The remaining routines were again scrutinized for syntax as related 
to a typical play of the game® From the observed functional and logical 
arrangement of these routines, the all important abstraction from *’real 
world,” SEALIFT, to model was attempted® 

Concurrent with, but not independent of, this abstraction, a variable 
which would be definitive to the expected value analysis was sought® Bear- 
ing in mind that the motive was to evaluate losses to the cargo fleet by 
submarine offensive, the natural and obvious choice was "the total expected 
number of torpedoes fired at cargo (vice convoy) ships during the war®" As 
stated, this number proved to be unattainable, the primary culprit being 
the time and event dependencies which permeate (are) the game® It appeared, 
however, that if the number of torpedoes a typical submarine would fire at 
cargo ships during a single cycle at sea, i®e®, from home port to station 
to home port, could be determined, this number could be extended into a 
total war figure® The intuitive appeal of this single cycle approach is 
that, save for the submarines initially on station, it removed, postponed 
really, the necessity for considering the element of time® It rendered 
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Submarine 
in port 



Submarine 

transits 

barrier 



Submarine 

transits 

barrier 




Submarine 
is on station 
XXXXXXXXXXXXXXXXX 




Convoy 




Legend 

► : Event flow 

XXXXXX : Convoy-submarine interaction 

FIGURE 2 

GRAPHICAL REPRESENTATION OF THE CYCLE ABSTRACTION 
FROM SEALIFT TO THE EXPECTED VALUE MODEL 
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individual submarines* unique and time dependent lives indistinguishable 
and it collapsed the immense and complicated time-event structure into a 
fairly tractable entity* The basic cycle is represented in Figure 2. The 
model for the expected value analysis developed in this thesis is therefore 
the logic needed, using the Markov chain process, to arrive at the number 
of torpedoes a submarine would fire at cargo ships during a single cycle 
at sea* 

In developing this logic, three difficulties became immediately apparent: 
1* Aside from the three different type submarines and the four pos- 
sible sonar boxes to which they could be assigned, the submarines initially 
on station are unique* These submarines have a probability distribution on 
their departure-f rom-on-station time, unlike other submarines which have a 
fixed mission time* This probability distribution necessitated a modifica- 
tion of the associated Markov chain, since prior SEALIFT results indicated 
this cycle to be a very sensitive factor* 

2. How are the results of a single cycle at sea "extended” to cover 
the total war? In particular, what assumptions are required to minimize 
the effects of a coarse approximation? How are excessively small, or large, 
convoy losses to be interpreted? 

3* How can the event and time dependencies within a single cycle be 
handled effectively? 

The third difficulty was resolved by essentially ignoring the event 
and time dependencies inherent to SEALIFT. For the most part, these depen- 
dencies proved incompatible with either the cycle abstraction or with the 
Markov chain, or both. Where possible, however, they were included within 
the framework provided by this thesis* 

The second difficulty was resolved by determining the minimum number of 
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cycles a live submarine would make during the ware The basic cycle was 
then repeated that minimum number of times with surviving submarines from 
one cycle serving as inputs to the following cycle* Rationale for this 
choice was provided by two facts* First, a priori information from the 
SEALIFT game indicated that where the submarine was almost immediately 
overwhelmed, the significant action took place early in the war* In this 
case, the number of cycles would prove immaterial* If, however, the con- 
voys were continually decimated, this expected value analysis would show 
abnormally high convoy losses, and this would actually serve the purpose 
intended by SEALIFT* 

If opposing forces were fairly balanced, successive cycles would show 
nearly identical convoy losses, again providing an indication of the progress 
of the war* In other words, under the assumptions that the omitted time and 
event dependent events are not significant, the analysis provided herein 
should reflect the trends of a corresponding SEALIFT game* This method 
of extension was also supported by the fact that, almost without exception, 
the omission of each of the dependent events effects an over-estimation of 
the number of cargo ship losses* The choice of a minimum number of cycles 
would then partially compensate for this over-estimation* 

The theoretical resolution of the first difficulty proved trivial, 
falling under the cognizance of the integral form of the theorem of total 
probability* The evaluation of this integral proved to be computationally 
messy* The functional simplicity of the probability distribution and a 
linear approximation to the step function integrand allowed a simple but 
accurate approximation to this integral* Each of these difficulties is dis- 
cussed more fully in subsequent subsections* 

Having chosen a tentative model to represent the game, this model was 
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divided into its component parts and separately analyzed. The basic .cycle 
could be broken into two phases, (1) the submarine-barrier penetration, 
and (2) the submarine-convoy interaction, the Markov chain of this develop- 
ment. The probabilities of the sequences of events associated with each 
phase, i.e., the barrier penetration and each state of the Markov chain, 
were analyzed via the corresponding 5EALIFT subroutines or parts of sub- 
routines. The SEALIFT subroutine logic was reformulated in the structure 
of a tree or truth table. While an attempt was made to construct the trees 
to match the subroutine logic exactly, some simplication was employed. 

Basic assumptions used in the expected value analysis . Perhaps the 
most definitive assumption made in the expected value analysis is that the 
submarine-convoy encounter should be the significant event for determining 
cargo ship losses. By implication, then, the convoy losses to air attack, 
both at sea and in port, and by barrier penetration are assumed small enough 
so as not to render the submarine target-limited, vice torpedo-limited, 
and also that each convoy be initially large enough to comply with this 
restriction. Thus, Event 3 (Convoy enters barrier CIC), Event 4 (Air attack 
against convoy). Event 6 (Convoy enters barrier C2C), Event 7 (Convoy in 
delivery port). Event 0 (Convoy enters barrier C2C on return trip), and 
Event 9 (Convoy enters barrier CIC on return trip) were classified as irre- 
levant to the analysis for this reason. 

This analysis hinges, then, on the number of submarines which penetrate 
the convoys, and prior to that, the number of submarines which reach the 
subline. For this reason, all the submarine loss events were deemed impor- 
tant. In the subsequent analysis, however, it was observed that the train 
of events which could occur to one submarine was dependent upon previous 
events involving other submarines. In other words, the stochastic process 
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associated with a particular submarine was dependent upon the realization 
of the process associated with other submarines on station with it and sub- 
marines which proceeded it to station. Such dependencies are inherently 
intractable to analysis and they were ignored in the expected value model. 

By their omission, these dependencies introduced a cumulative over-estimation 
of the number of cargo ship losses. In particular, Event 17 (Submarine enters 
variable barrier) calls the subroutine VARDET to recalculate the reduced 
barrier-detects-submarine probability in the event a submarine sinks one 
or more of the barrier ships. Following submarines would then stand a better 
chance of slipping through that barrier until the loss was replaced (Event 22 
(Replace e barrier ship)). A concurrent effect is the torpedo usage on 
barrier ships. This expenditure directly affects the ammunition available 
to be used on cargo ships later in the submarine’s cycle. This last effect 
could become more invidious during a HUK attack. Event 15 (Submarine versus 
HUK group). In this routine the only ways a submarine could escape, if de- 
tected by the HUK Group, are if her time or ammunition runs out or if her 
torpedoes miss their mark! In addition to introducing event dependencies, 
the HUK attack event presented further theoretical difficulties in the form 
of time dependencies. Each submarine, while at sea, is subject to HUK attack 

t 

every HUKHR interval. This event meshes with every other event in the sub- 
marine’s life and tends to hopelessly confuse expected value analysis. 

The subroutines FSEDET and ASWDET , which perform the same function for 
the forward screen escort and the ASW screen as VARDET does for* the barrier 
ships, are called during the play of Event 19 (Submarine versus convoy). 

Like VARDET, these routines influence the train of events of following sub- 
marines. Again, parasitic torpedo usage is encountered; however, an attempt 
was made to approximate this usage in the model in the case of the forward 



15 







& 

I 






screen escort 



Another, and time related, dependency was the possible monthly change 
in submarine-detects-convoy probability for the submarines stationed in the 
fourth sonar box, a portion of Event 5 (Convoy arrives at submarine line). 
This implied an incompatibility with the cycle concept. Only fortuitous 
combinations of the number of cycles of the war, duration of the war and 
constancy of the detection probability associated with the fourth sonar box 
during each cycle would permit easy inclusion of this dependency. 

Similarly, Event 13 (Strike at submarine source) was a time dependent 
event. Provision for up to ten air strikes, at times chosen through the 
input parameter set, were made in SEALIFT. These air strikes, too, are 
incompatible with the cycle concept used in the model. 

To summarize, the routines ASWDET , F5DET, VARDET, Event 15 (Submarine 
versus HUK group), Event 22 (Replace a barrier ship). Event 13 (Strike at 
submarine source), and a portion of Event 5 (Convoy arrives at the submarine 
line) were ignored by virtue of their intractability to analysis. 

Finally, Event 12 (Build submarines) provides for a monthly submarine 
b&ilding rate. While this event is also time dependent, it was adapted to 
the model by applying the number of submarines built during the entire war 
as an additional input to the first cycle, a striking illustration of the 
time independence of the cycle concept. 

It should be noted that, in SEALIFT, there exist at least two levels of 
logic, the internal event logic which computes losses to convoys, submarine 
losses, etc., and the logic which controls the flow of the game, i.e., com- 
putes when and how events will occur. These levels are frequently meshed 
within each subprogram, and it is with regard to the first named level that 
routines were ignored or classified as irrelevent. The second type logic 
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was used, implicitly, whenever necessary, in the analysis® In addition, 
the following were ignored as being primarily bookkeeping type routines: 
Event 1 (Build ships), Event 2 (Form convoy). Event 10 (Convoy dispersed 
at home). Event 14 (Submarine leaves port). Event 10 (Submarine arrives on 
station). Event 20 (Submarine goes home). Event 21 (Submarine arrives home), 
and Event 24 (Replace convoy ships)* 

The events from SEALIFT that form the basis for the analysis are there- 
fore: Event 5 (Convoy arrives at submarine line). Event 11 (Initialize 

submarines). Event 16 (Submarine enters constant barrier). Event 17 (Sub- 
marine enters variable barrier), and Event 19 (Submarine versus convoyK 

3*1 Markov Chain Structure for the Expected Value Analysis 

For a given submarine, the initial Markov chain begins after that sub- 
marim^s first convoy has passed the subline. The initial states for this 
Markov chain are the possible sequences of events which could have occurred 
to the submarine when that convoy passed. These states may be represented 
by the arborescence , with root ’’Submarine on Station,” exhibited in Figure 3 
where a state of the Markov chain is defined as the path from the root to a 
terminal point. Letting L^ denote the states of interest, define 

E^ = submarine does not detect the convoy. 

E 2 = submarine detects the convoy, does not penetrate the Forward 
Screen Escort (F5E), and is not killed* 

E^ = submarine detects the convoy, penetrates the F5E, penetrates the 
inner screen (A5W), shoots one salvo at cargo ships, and survives 
any resulting Search and Attack Unit (5AU) action* 

E^ = submarine detects the convoy, penetrates the F5E, penetrates the 
A5W, shoots one salvo at cargo ships, and is killed in resulting 
5AU action. 
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FIGURE 3 

ARB0RE5CENCE ASSOCIATED WITH THE STATES OF THE FIRST MARKOV CHAIN 
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= submarine detects the convoy, penetrates the FSE, penetrates the 
A5W, shoots two salvos at cargo ships* and survives any resulting 
SAU action* 

= submarine detects the convoy, penetrates the FSE, penetrates the 
ASW, shoots two salvos at cargo ships, and is killed in the re- 
sulting SAU action* 

= submarine detects the convoy, penetrates the FSE, and is killed 



by the ASW. 



E 
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= submarine detects the convoy and is killed by the FSE* 



This Markov chain allows for the possibility that if a single salvo 
does not exhaust the torpedo supply of a submarine during an attack, there 
is a chance that the submarine will be unable to fire again at that convoy 
and will, therefore, return to station with a partial load of torpedoes to 
await the arrival of the convoys* It is assumed in this model, however, 
that two salvos will exhaust a submarine’s torpedo supply* This assump- 
tion is equivalent to the assumption that the necessary train of events which 
would place a submarine in the position of having fired more than two single 
salvos at as many different convoys to be very improbable, and neglects this 
possibility* ^ 

The occurrence of states E^^ and E^ imply that the submarine returns to 
station where she awaits the next convoy arrival, which is to say, she re- 
mains alive in this Markov chain* State E^ is followed by the submarine’s 
attempt to return home, and is, therefore, an absorbing state* Similarly^ 

states E., E^, E-,, and E„ are absorbing since the submarine is killed* State 
4 b I o 

E^ is absorbing for another reason; it is the entry point for the second 
Markov chain described below* 

If the probability of the initial states are denoted by Pr(E^) = p^, 
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the one step transition probability matrix is given by equation (2*4) . 

If a total of n+1 convoys pass the subline during a submarine’s allotted 
on-station time, then the transition matrix, P, must be raised to the n th 



power since the initial states are caused by passage of the first convoy to 

yield the n step transition matrix The resulting matrix, P^, is given 

by equation (2*6). This matrix and the initial probability distribution, 

p^, suffice to determine the unconditional probabilities, that a 

submarine will be in each of the possible states at the completion of n 

steps, (n+1) convoys having passed* The unconditional probabilities, 

are obtained from equation (2*3) and from them, the probabilities of certain 

relevant events can be obtained: 
d 



(3^0) 



P(A^) 



= Pr (a submarine survives this Markov chain to attempt 

passage home/ n+1 convoys have passed the subline) 

p(n) (n) (n) 

1 2 5 



(3.1) 



(3.2) 






PCAj) 



Pr (a submarine fires one salvo at the cargo ships in 
this Markov chain/ n+1 convoys have passed the 

subline ) 

p(n) p(n) 

3 4 

Pr (a submarine fires two salvos at the cargo ships in 

this Markov chain/ n+1 convoys have passed the 

subline) 

. p'"' 

b b 



(3.3) 



P(B ) 
r 



= Pr (a submarine is absorbed in state upon passage 

of the r th convoy) 
r-1 



P3 (PjL + P2> 



r “ 1....^ n ■4' 1 
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The second Markov chain is almost identical to the first chain and 



is concerned with those submarines which are absorbed in state E^o These 
submarines return to station having fired a single salvo of torpedoes at 
the previous convoy « Assuming that a total of (n*fl) convoy cross the^ sub- 
line, the number of steps of the second Markov chain depends upon which 
convoy caused a submarine to be absorbed in state E^* If a submarijie is 
absorbed in state E^ at the passage of the r th convoy, the {r+l)st convoy 
"initializes" the second Markov chain and there are (n*»-l) - (r+2) 4- 1 = n-r 
steps in the second chain® The probability that a submarine enters the 
second Markov chain "just after" passage of the r th convoy is given by 
equation (3«3)o 

The second Markov chain is described by its states: 

= submarine does not detect the convoy® 

= submarine detects the convoy, does not penetrate the F5E , 
and is not killed® 

F^ = submarine detects the convoy, penetrates the FSE , pene- 
trates the A5W, shoots her remaining torpedoes at cargo 
ships, and survives any resulting SAD action® 

F^ = submarine detects the convoy, penetrates the FSE, pene- 
trates the A5W, shoots her remaining torpedoes at cargo 
ships, and is killed in the resulting SAU action® 

F^ = submarine detects the convoy, penetrates the FSE, and is 
killed by the ASW® 

Fg = submarine detects the convoy, and is killed by the FSE® 

Denoting the probabilities of the initial states of this Markov chain by 
Pr(F.) = q., the one step transition probability matrix is then 



21 



II 

1 



I 



Q = 



‘’i ‘’z ’^3 ‘’4 ‘’s 



qi ^2 '’a '’4 '’s '’e 

001000 

000100 

0 0 0 0 1 0 

000001 



= (q . . ) 
ij 



and the n-r step transition matrix, by analogy with equation (2,6), is 



r.n-r 


, .n-r -1 

+ q2> 


qi q2‘ 


, t \n-r 

1 - (qi + qz^ 


P 3 P 4 P 5 ^ 6 ' 


Q 




. X 








qi q2 


1 - (qi + q2) 


P 3 P 4 P 5 P6_ 



, (n-r). 
‘■’ij > 



4x2 



4x4 



,(n) 



If is defined as the unconditional probability that the submarine is in 

state F. after n steps, then, using the theorem of total probability, 

J 



J 



n+1 

z 

r=l 



n 

z 

r=l 



Pr (a submarine is in state F. after n+1 convoys/ 

J 

she entered the second chain after the r th 
convoy has passed) x Pr (she entered the second 
chain after the r th convoy). 



X P(B^) , n ^1 



^ ^ (n-r) , .r-1 > , 

.1 ‘’i ‘^ij P3 ‘Pl ^ ^ 



r=l i=l 



where 



(o) d 

'ij 



6ij. 



the Kronecker delta. 
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The probabilities of the remaining events of relevance can now be ob- 



tained as: 
(3o4) 



and 

(3.5) 



P(A^) = Pr (a submarine survives the submarine-convoy inter- 

action/ n+1 convoys have passed the subline )< 

¥ 

(n) (n) (n) (n) „(n) (n) . ^ r. 

= P^ + P^ + P5 + + R^ + R3 + P3(Pi + P2> 



P(Ag) = Pr (a submarine fires her remaining salvo in the second 
Markov chain/ n+1 convoys have passed the subline) 






In SEALIFT, the submarines may be separated into two sets, those which 
are on station at the beginning of the war and have a departure time distri- 
buted uniformly on an interval, and those submarines which arrive on station 
after the war starts and have a fixed mission timeo Once the game para- 
meters have been set the arrival of each convoy is deterministic® As a 
consequence, the variable n+1 is a step function of the submarines departure 
time T, ioe®, n+1 = g(T)® Representing the equations, (3oD), (3«1), (3®2), 
(3o4) and (3® 5) symbolically as 

Pr(A/n+l) = Pr (of the event A given n+1 convoys have passed) 
this probability becomes Pr(A/g(T) )® Invoking the theorem of total probabil- 
ity, Pr(A) can be found as 



(3®6) 



Pr(A) 



/ 



Pr(A/g(t)) d Pr(T^t) 



This equation is valid for both sets of submarines although trivial for 
those with a fixed mission time® For the submarine with a fixed mission time 
Pr(A) is simply the integrand of this equation evaluated at the upper limit 
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of to The evaluation of equation (3o6) is crucial in determining the ex- 
pected number of cargo ships sunk per cycle. 

This completes the construction of the Markov chain for a generalized 
submarinBo To apply this structure to the SEALIFT game, the integrals Pr(A) 
and the transition probabilities p^^ and must be expressed in terms of 
of the FORTRAN variables of SEALIFT, The submarine barrier penetration must 
be analyzed in like terms* Finally, the evaluated terms must be synthesized 
to answer the basic question ”What is the expected number of torpedoes fired 
at cargo ships per cycle?” and its extension ”What is the expected number of 
torpedoes fired at cargo ships in the war?” 

3o2 Detailed Analysis in Terms of the FORTRAN Variables of the SEALIFT Game* 
The detailed analysis to follow will be executed in the following order, 
(1) evaluation of the integrals Pr(A), (2) evaluation of the barrier penetra- 
tion, the transition probabilities, and related expected number, (3) syn- 
thesis to determine the expected number of torpedoes fired at cargo ships 
during the war. The following list of SEALIFT I ^ij FORTRAN variables is 
pertinent to the analysis, and the following convention is observed: The 

letter ”1” is reserved for the type submarine index; however, it is employed 
in two ways* As a subscript. Bog,, ALTIM (I), ”1” can take on the values 1, 
2, and 3. As a part of a variable name, e.g*, NSUB*I*T, ”1” is set off in 
this analysis by single quotation marks and can take on the values 1, 2 and 
N* In both usages, the values of ”1” denote conventional number 1, conven- 
tional number 2, and nuclear submarines, respectively* 

Variable Name Definition 

TBCON Time between each convoy sailing (hours)* 

CONSPD Speed of convoy (knots)* 
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CONHR 


Time after start of war that the first convoy leaves 
home port in hours® 


ONEWAY 


One way convoy distance from home port to delivery port 
(n® miles)® 


ULTIME 


Time convoy spends in delivery port. 


ALTIM(I) 


Allowable time in hours for type I submarine before it 
must leave its station® 


DSUBTA 


Distance from submarine home port to submarine station 
(n® miles)® 


SUBSPD(I) 


Speed of type I submarine® 


NSUB'I'T 


Number of type I submarines at start of war® 


FSUB'I'OS 


Fraction of type I submarine on station at start of war® 


P'I'BOX(J) 


Probability that a type I submarine on station will be 
in BOX # J, J = 1, 4. 


NSUB'r 


Number of type I submarines built per month. 


PBARSCd ,J) 


Probability that a type I submarine is sunk by the Jth 
constant barrier, J = 1, ®.®, 4. 


PBAR5V(I,J) 


Probability that a type I submarine is sunk by the Jth 
variable barrier, given it was detected by this barrier, 
J “ 1, ••., 4® 


PVBDET(I,J) 


Probability that a type I submarine is detected by the 
Jth variable barrier, with barrier at its full quota of 
ships® 


VAIRDET(J) 


Part of detection probability of variable barrier J 
that is not supplied by ships within that barrier which 
might be sunk by a submarine® Thus, if a mixed barrier 
of aircraft and ships were used and all the ships were 
sunk, the probability of detection drived from the air 
barrier (VAIRDET ( J) ) would still be present. The ”total" 
detection probability of a variable barrier is PVBDET(I,J) 
+ VAIRDET (J)® For a ship barrier only, VAIRDET(J) should 
be set to zero® 


RGDT(K,I) 


Maximum athwartship detection range of type I submarine 
against convoy type K sonar conditions 
K=1 iff good sonar conditions 
K=2 iff poor sonar conditions 


SUBLINE 


Length of submarine line (n® miles), i.e., width of 
patrol area® 


TWAR 


Duration of war (hours). 
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FAIRSC 


Is forward escort screen made up entirely of aircraft 
(non-sinkable )? (YES=1, N0=2)« 


ATSFE(I) 


Number of times a type I submarine can try to repenetrate 
forward screen 


PDFS(I,K) 


Probability that a type I submarine in type K sonar con- 
ditions (denoted by **(I,K) submarine” henceforward) is 
detected by forward screen escort while penetrating 
screen with screen at its full quota of shipso 


PDFER 


Part of detection probability of forward screen that is 
not supplied by ships within that screen which might be 
sunk by a submarine* Thus, if a mixed screen of air- 
craft and ships were used and all the ships were sunk, 
the probability of detection derived from the air screen 
(PDFER) would still be present. The “total” detection 
probability of the forward screen is then PDF5(I,K) + 
PDFER. For a ship's screen only, PDFER must be set to 
zero. 


PDFCE(I,K) 


Probability that an (I,K) submarine is sunk by forward 
screen escort, given submarine has been detected by same. 


PSF5E ( I ) 


Probability that a type I submarine shoots at forward 
screen escort ship in an attempt to break contact. 


PFCES(I,K) 


Probability that an (I,K) submarine sinks forward screen 
escort ship, given it decides to shoot at it. 


PSFCEd ,K) 


Probability that an (I,K) submarine is sunk by forward 
screen escort ship, given this submarine has already sunk 
a forward screen escort ship, i.e. , flaming datum. 


ASALVO(I) 


Number of anti-ASW torpedoes used per salvo by a type I 
submarine. 


PDETUd.K) 


Probability that an (I,K) submarine is detected by the 
ASW inner screen with direct path sonar, given it was 
not detected previously with bottom bounce sonar. 


VA(K) 


Conversion factor reflecting improved ASW inner screen 
detection probability with direct path sonar, in type K 
sonar conditions, because of previous bottom bounce de- 
tection. 


PBB(I,K) 


Probability that an (I,K) submarine is detected by the 
ASW inner screen with bottom bounce sonar. 


PKBB(I,K) 


Probability that an (I,K) submarine is killed, given 
that it is attacked or re-attacked following a detection 
with bottom bounce sonar. 
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PAT(I,K) 

PCC(I,i) 

PABBd ,K) 
PKBBG(I) 
QTC 177 

PKINS(I) 

PSCA(I,i) 

PAA(I) 

PKSC(I) 
PSUBAT(l) 
PKOSC(I) 
QTU 390 



Probability that an (I,K) submarine is attacked given 
that it has been detected with bottom bounce sonar only. 

Probability that a group from the ASW inner screen re- 
attacks a type I submarine at the i th opportunity in 
action resulting from a bottom bounce detection, 
i=l, 2, 3o 

Probability that an (I,K) submarine attacks a group 
from the ASW inner screen in action resulting from a 
bottom bounce detection. 

Probability that a type I submarine sinks a ship in 
the group from the ASW inner screen in action resulting 
from a bottom bounce detection. 

Designator for possibility that an engagement between 
the ASW inner screen and a submarine is terminated after 
an unsuccessful submarine counter attack, e.g., due to 
quiet torpedoes used by submarine (QTC equal to zero 
implies engagement is to continue; QTC equal to one 
implies engagement is terminated. No other values are 
possible) o 

Probability that a type I submarine is sunk by the ASW 
inner screen during action following a detection with 
direct path sonar. 

Probability that an ASW inner screen re-attacks a 
type I submarine at the i th opportunity in action re- 
sulting from a detection of sub by ship with direct 
path sonar, i=l, 2, 3. 

Probability that a type I submarine attacks the ASW 
inner screen after an unsuccessful attack on the sub- 
marine by the screen resulting from a detection with 
direct path sonar. 

Probability that a type I submarine kills an ASW inner 
screen ship in action resulting from detection of sub 
by ship with direct path sonar. 

Probability that a type I submarine attacks the ASW 
inner screen following an unsuccessful bottom bounce 
detection attempt by the screen. 

Probability that a type I submarine sinks an ASW inner 
screen ship without prior detection of the submarine by 
the screen. 

Designator for radiated noise level of enemy anti-screen 
ship torpedoes. (QTU equal to one implies a quiet tor- 
pedo; QTU equal to zero implies a noisy torpedo. No 
other values are possible). 
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VM{K) Conversion factor reflecting change in the ASW inner 

screen detection probability with direct path sonar in 
type K sonar conditions because of flaming datum result- 
ing from sinking of screen ship by previously undetected 
submarine o 

CSALVDd) Number of anti-shipping torpedoes used by type I sub- 

marine« 



PD5AU(I,K) 



PSSAU{I,K) 

ATTSAU(I) 



Probability that an (I,K) submarine is detected by the 
5AV, given the submarine has sunk one or more cargo ships 
in the convoy® 

Probability that an (I,K) submarine is sunk by the SAL), 
given the SAD detected the submarine* 

Does a type I submarine attack a number of the SAU group, 
given it was not sunk by this group, YES=1* 



RASW(I) 



Number of times a type I submarine can re-attack the 
convoy after SAU action* 



TORPS(I,ASW) Number of anti-ASW torpedoes per type I submarine* 

Note: If only one type of torpedo is used for the play 
of the game and it is desired to fire this torpedo at 
all targets, then it should | must 1 be stored in TORPS 
(I, ASW). L J 



3*2*1 Evaluation of the Integrals Pr{A) 

To evaluate the integrals Prt"A)^> the equivalences between the first and 

second Markov chains are exploited to condense the notation. The states E^^, 

E_, E_, and Eq are seen to be equivalent to the states F, , F_, F-, and F^ 

2 f o 12 3 0 

respectively, and therefore, p. = q. , p_ = q_, p_ = q_ and p. = p.* For 

1122f3 bo 

3 = (pl + p2^ 



and 



1-s 



n+1 



f (n) = 



1-s 



and letting D denote the derivative operator with respect to 
tion of the unconditioned probabilities, arid yield 
relationships: 



s, the evalua- 
the following 
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ii 



,(n) 

1 


B 


Pl» 




,(n) 

2 


= 


P2= 




,(n) 

j 


= 


Pjf(n), j = 




= 


Pi P3 


r. n 
Os 






P2 P3 


Df(n) 


,(n) 

‘3 


= 


P3 P3 


Df(n) 


,{n) 

‘4 


= 


P4 P3 


Df(n) 


,{n) 

5 


= 


P7 P3 


Df (n) 


(n) 

6 


= 


P8 P3 


Df(n) 



Using these values, equations (3*0, (3ol), (3o2), (3o4) and (3, 
section 3ol are evaluated to yield 



(3.0)' 


P(Aj^) 


n+1 

= s 


^ P 5 


f (n) 


(3.1)' 


P(A2> 


= ^P3 


+ P 4 J 


f (n) 


(3.2)' 


P(A3) 


= ^P5 


+ Pg) 


f (n) 


(3.4)' 


P(A^) 


= (1 + 


P 3 D) 


n +1 

8 


(3.5)' 


P(A^) 


= ^^3 


^ P4^ 


P 3 D 



Generically each of the integrals Pr(A) 
Pr(A) = J Pr(A/g(t)) d Pr(T 



= t) 



is a function of the integral N 

,g(t) 



N = 



L 



d Pr(T t) 



5 ) from 
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and the first derivative of N with respect to s, assuming the integration and 

differentiation operations can be interchanged. The distribution Pr(T = t) 

is non-trivial only for the submarine initially on station and is, in SEALIFT, 

dt 



uniform on the interval (0, ALTIM) , i.e., dPr(T t) = 



ALTIM 



To explore the structure of the function g(t) and remain within the 
FORTRAN concept of SEALIFT, two new FORTRAN variables are defined 

TCON(I) = The game at which the I th outgoing convoy crosses 
the subline, 



and 

RTCON(I) = The game time at which I th returning convoy crosses 
the subline. 

For those submarines initially on station, g(t) is zero for t < TCON(l), and 
steps to one at t = TCON(l). It then increases by one after each interval 
TBCON until the first returning convoy crosses the subline, at game time 
RTCON(l). Until the SEALIFT game parameters have been set, the fine structure 
of g(t) is not well defined after RTCON(l). It does, however, increase by 
two every TBCON hours after RTCDN(l). Clearly then, the integrals Pr(A) 
involve extensive summations. To avoid these summations a piecewise contin- 
uous linear approximation, h(t), to the step function, g(t), is sought such 
that for g(t) = h{t), the value of N is easily found, viz.. 



(3.7) 



N 




dt 



h(t) 

a 3 

Ln(s) dh( t ) 
dt 



Incorporating the SEALIFT FORTRAN variables and the fact that in SEALIFT 
the subline is positioned half way between convoy home port and delivery port 
it can be seen that 
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TCON(l) 



+ COIMHR 



ONEWAY 
2 X CONSPD 

T RNTWAY 

RTCON(l) = Y X + CONHR + ULTIME 

and further, 

TCON(I) = TCON(l) + (I-l) TBCON 
RTCON(I) = RTCON(l) + (I-l) TBCON 

If K is defined as the number of outgoing convoys which pass the subline 
before the first returning convoy passes the subline, then K is the largest 
integer which satisfies the inequality 
TCON(K) 5 RTCON(l) 
or 

(K - 1) TBCON + TCON(l) ^ RTCON(l) 

+ 1 , 

where the brackets denote "the greatest integer less than or equal to." 
Introducing a new FORTRAN variable H(I) for the expected value analysis, the 
values of H(I) and the approximation h(t) are defined as: 



H(l) 


d 


h(0) = 0 






H(2) 


d 


h(TC0N(l) - 


0) = 0 




H(3) 


d 


h(TC0N(l) + 


II 

o 




H(4) 


d 


h(RTCONd) • 


■ ^ TBCON 


(RTCON(l) - TCON(D) 


H(5) 


d 


h(RTCONd) ■ 


¥ 0 ) = K + 




H(6) 


d 


h(ALTIM) = 


(K + t) + 


(ALTIM - RTCON(D) 



Figure 4 depicts both a sample step function g(t), and its linear approxi- 
mation h(t). Using h(t), the generic integral, N, becomes 



RTCON(l) - TCON(l) 
TBCON 
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(3.8) 



N = 



'TCON<D /RTCON(l) /ALTIM 

s dt + / s dt + / 8 dt 



/TCON(l) 



(RTCON(l) 



and, by equation (3,7), becomes 



H(4) H(3) 

= TCON(l) 4- TBCON . 

Ln ( s ) 2 



H(6) H(5) 

a =_S 

Ln(s) 



The following derivatives are of immediate interest in evaluating the inte- 
grals Pr(A)^ 



^ N 
ALTIM 



xnrnM . ( \ Iu(a\ H(4)-1 H(3)-l. , H(4)-l H(3)-l. 

TBCON Ln(s) (H(4) 8 - H(3) s ) - (s -s ) 



ALTIM 



Ln^(s) 



rnroM ^ H(6)-l H(5)-l, , H(6)-l H(5)-l, 

TBCON Ln(s)(H(6)s -H(5)s )-(s -s ) 



2 ALTIM 



Ln^(s) 



and 



ALTIM - N 
ALTIM (1 - s)^ 



1 - 



N 



ALTIM 



- (1 - s) D j 



^1 



ALTIM 



(1 - s)‘ 



In terms of these values, the integrals Pr(A) become 

(3*9) PAl i Pr (a submarine initially on station survives the first 

Markov chain to attempt passage home) 



N ALTIM - N 

ALTIM ^5 1 - s 

(3«10) PA2 = Pr (a submarine initially on station fires one salvo at 

the cargo ships in the first Markov chain) 



, . ALTIM - N 

'^3 * ^4^ ALTIM (1 - s) 
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FIGURE 4 

A SAMPLE STEP FUNCTION, g(t), AND ITS APPROXIMATION, h(t) 
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(3.11) 



PA3 = Pr (a submarine initially on station fires two salvos at 
the cargo ships in the first Markov chain) 



, , ALTIM - N 

^^5 ^6^ ALTIM (1 - s) 



(3ol2) PA4 = Pr (a submarine initially on station survives the inter- 
action) 



^3 ALTIM 



+ (Pg + Pg qg D) 



ALTIM - N 
ALTIMd - s) 



(3,13) PAS = Pr (a submarine initially on station fires on cargo ships 

in the second Markov chain) 



= (qg + q^) P3 D 



ALTIM - N 
ALTIM (1 - s) 



This completes the evaluation of the integrals Pr(A) for those submarines 
initially on station® 

The integrals Pr(A) are trivially Stieltjes integrals for those sub- 
marines which arrive on station after the war starts. The following is based 
on the assumption that the arrival time for these submarines is later than 
RTCON(l), ioBo, the full convoy arrival rate of two convoys per TBCON hours 
has been achieved prior to submarine arrival at the subline. In this case, 



Pr(A) = Pr(A/ g(maximum value of t)) 



but 



g(max t) = 



ALTIM - submarine transit time to station 
TBCON 



X 2 



Neglecting possible delays engendered during HUK attacks and during the barrier 

D5UB5TA 

penetration, the submarine transit time is S ~ ubSPd" ' therefore 
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1 



Pr(A) = Pr(A/ n+1 = 



ALTIM 



DSUB5TA 

SUB5PD 



TBCON 



X 2 



Letting KONE = 



be found 
(3.14) 



ALTIM - 



DSUBSTA 

SUBSPD 



TBCON 



X 2 



, the following probabilities can 



PA 6 = Pr (a submarine not initially on station survives the first 
Markov chain to attempt passage home) 

KONE 



KONE 1 - s’ 

= j-TT 



(3,15) 



(3.16) 



(3,17) 



(3.10) 



PA7 = 



Pr (a submarine not initially on station fires one salvo 
at the cargo ships in the first Markov chain) 

KDNE 



(P3 + P4) 



1-3 

1-5 



PA 8 = 



d 

PA9 = 



PAIO = 



Pr (a submarine not initially on station fires two salvos 
at the cargo ships in the first Markov chain) 

, KONE 

/ \ 1-3 

•Ps ^ P6’ TTT 

Pr (a submarine not initially on station survives the sub- 
marine-convoy interactiohj 

vnMr 1 KONE 

l^ r>\ KONE , 1 - s 

(1 + P 3 D) s + (Pg + P 3 P 3 D) 

Pr (a submarine not initially on station fires on cargo 
ships in the second Markov chain) 

KONE 

‘■’3 * P.’ Pa “ rH 



This completes the evaluation of the integrals in terms of the evaluated or 
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known FORTRAN variables and the as yet unevaluated transition probabilities. 

3o2o2 The Barrier Penetration, Transition Probabilities and Related Expected 
Numbers 

Until now, it has been convenient to suppress in the analysis any nota- 
tion concerning submarine type and sonar condition. The Markov chain notation 
was difficult enough without the burden of two additional subscripts. The 
remainder of this analysis will distinguish submarines by type and by assigned 
sonar condition. It will be seen that each of the probabilities given by 
equations (3.9) through (3.18) is a function of submarine type and sonar 
condition. 

Placement of submarines . In 5EALIFT, there are eight ways to classify 
submarines of a particular type by placement at the start of war: two choices 

for on station or not and four choices of sonar box. This classification is 
done by Monte-Carlo methods. Type sonar condition within each box is de- 
termined by input. For the purposes of this analysis, sonar box classifica- 
tion is immaterial except for the restriction that the fourth sonar box is 
not included in the analysis. For the initial placement of submarines the 
expected number of type I submarines initially on station in type K sonar 
conditions is given by 

E2(I,K) = NSUB'I*T X FSUB*I*0S X 

where the summation is over all sonar boxes assigned sonar condition K. The 
expected number of type I submarines not initially on station but assigned 
to type K sonar conditions, E3(I,K), is derived from two sources, the fraction 
of on hand submarines not initially on station and the submarines built during 
the war. Using a nominal thirty day month, 

E3(I,K) = (NSUB’I’T (1 - FSUB’I'OS) +NSUB’!' x x 5]p’I'B0X(J) 
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where the summation is over all sonar boxes assigned sonar condition K® 



Barrier penetration ® A submarine will successfully transit the eight 
distinct barriers if and only if she successfully transits each barrier. In 
the SEALIFT game, the barriers are evenly divided between two types, constant 
and variable® The latter type are variable in the sense that ^the barrier- 
detects-submarine probability will vary with barrier strength® Under the 
basic assumptions used in this expected value analysis, however, each of 
the barriers is treated as a constant barrier® For constant barriers there- 
fore 

PA11(I,J) = Pr (a type I submarine successfully transits the J th 

constant barrier) 

=5 1 - Pr (she is sunk in this barrier) 

= 1 - PBAR5C(I,J) 



and 

PA12(I,J) = Pr (a type I submarine successfully transits the J th 

variable barrier) 

= 1 - Pr (she is sink in this barrier) 

= 1 - PBARSV(I,J) X (PVBDET(I,J) + VAIRDET(J)) 



Therefore the probability that a type I submarine successfully transits the 
barriers, P0(I), is 



4 

P0(I) = 7T PA11(I,J) X PA12(I,J) 

J=1 



Convoy detection ® In the SEALIFT game, all convoys cross the subline 
at the subline’s midpoint® Detection probabilities are of the "cookie 
cutter" type where submarines are placed (at each convoy crossing) uniformly 
on the subline® The probability an (I,K) submarine detects a convoy on any 



36 



I 



one trial, 1 - P1(I,K), is then 



1 - PX(I,K) 



2 X RGDT(K,I) 
SUBLINE 



The transition probabilities, " P1(I,K), are thereby fully defined, 

and are dependent on submarine type and on sonar condition* 

Penetration of the forward screen escort * This event and the following 
events are the substance of Event 39 (Submarine versus convoy) of the SEALIFT 
game* Figure 5 is the simplified truth table (tree) of that part of Event 19 
from which the F5E associated probabilities and expectations used in this 
analysis are derived* All of these probabilities are conditioned on the 
submarine having detected a convoy* The input parameters FAIRSC and ATSFE(I) 
are obviously not probabilities* FAIRSC is transformed into a 0,1 parameter 
in this analysis by the transformation 2 - FAIRSC, allowing its use in prob- 
ability statements* ATSFE(I) is another tactical limitation of the type I 
submarine and this parameter is usually a "small** integer for the conventional 
submarine and a **large** integer for the nuclear submarine. If ATSFE(I) = 0, 
a type I submarine will not be allowed to attempt repenetration of the FSE 
in the event she is rebuffed* For this submarine in type K sonar conditions 
the following probabilities are obtained 



P7(I,K) 



d 



Pr (a submarine penetrates the FSE; ATSFE(I) = 0) 



= 1 - (PDFS(I,K) + PDFER) 

P9(I,K) = Pr (a submarine returns to station; ATSFE(I) = 0) 



(PDFS(I,K) + PDFER) x (1 - PDFCE(I,K)) x (1 - (2 - FAIRSC) 



X PSFSE(I) X PFCES(I,K) x PSFCE(I,K)) 



and 



P6(I,K) 



d 



Pr (a submarine is killed attempting FSE penetration; ATSFE(I)«0) 



1 - P9(I,K) - P7(I,K) 
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penetrates F5E 
returns to station 
is killed 



FIGURE 5 

ARBORESCENCE ILLUSTRATING THE SUBMARINE-FORWARD SCREEN ESCORT INTERACTION 
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In general, if ATSFE(I) = n, n attempts at penetration will be allowed, and 
P7(I,K,n) = Pr (a submarine penetrates the FSE; ATSFE(I) = n) 



= P7(I,K) 



1 - (P9(I,K)) 
1 - P9(I,K) 



n+1 



and 



P9(I,K,n) = Pr (a submarine returns to station; AT5FE(I) = n) 

= (P9(I,K))"^^ 



P8(I,K,n) = Pr (a submarine is killed attempting FSE penetration; 

ATSFE(I) = n) 



_ pR(T K) 1 - (P9(I>><)) 
- 1 - P9(1,K) 



n+1 



As derived herein these values are predicated only on the approximation that 
the F5E-detects-submarine probability has not been degraded by prior sub- 
marine attack either by the submarine under consideration (ATSFE{I)> 0) or 
by earlier submarines and that when the submarine attempts to attack the 
screen she will have sufficient torpedoes to do so« To find the expected 
number of torpedoes a submarine fires at FSE ships given she penetrates the 
screen, suppose ATSFE(I) = n, and define the following events 

Cj = event that a submarine penetrates the screen on the j th 
try, j = 1, o o • , n^-1 

= event that a submarine fires i salvos, each of size A5ALV0(I) 
during an attempted penetration, i = □, •••, n+1 
Then El(I,K) can be defined as the expected number of torpedoes fired by an 
(I,K) submarine during an attempted penetration of the FSE given she penetrates 
the FSE and ATSFE(I) = n, ioeo, 
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(3.19) 



n+1 

n Z! 

E1(I,K) = ASALVOd) i 1 

i=0 P7(I,K,n) 



The probability Pr (C^, D^) can be derived by the following argument: the 

event ^ sequence of j trials where the j th and last trial is 

a non-detection by the FSE® This trial is preceeded by i trials (not 
necessarily in sequence) in which the submarine is detected, not sunk, does 
attack the F5E , and either sinks an F5E ship and is not killed or does not 
sink an F5E ship. The last trial is also preceeded by j - i - 1 trials in 
which the submarine is detected, not sunk, and does not attack the F5E. 

Denote the probability of each type trial, in the order listed, by Tl, T2, 
and T3. Then 

T1(I,K) = P7(I,K) 

T2(I,K) = (PDF5(I,K) + PDFER) x (1 - PDFCE(I,K)) x (2 - FAIR5C) 

X P5F5E(I) X (1 - PFCE5(I,K) x PSFCE(I,K)) 

and 

T3(I,K) = (PDF5(I,K) + PDFER) x (1 - PDFCE(I,K)) x (1 - (2 - FAIR5C) 

X PSFSEd)) 

Therefore , 



Pr(C., D.) = 

‘ J 1 



j-1 

i 



X (T2(I,K))^ X (T3(I,K))-^"^”^ x P7(I,K) 



and equation (3.19) becomes 



(3.19)' 



E1(I,K) 



ASALVO(I) 




i=l 



n+1 




(T2(I,K))^ X (T3(I,K))-^"^"^ 

1 - (P9(I,K))' ^'*'^ 

1 - P9(I,K) 
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Penetration of the inner screen (ASW)* Having penetrated the FSE , an 
attacking submarine must penetrate the inner screen prior to reaching the 
cargo shipso Figure 6 is the simplified tree of that part of Event 19 
(Submarine versus convoy) of SEALIFT from which the ASW associated probabil- 
ity of penetration is derived for this expected value analysis* This portion 
of the analysis is based upon an undated and informal modification of Event 
19 which, however, modifies only the ASW penetration portion* The assump- 
tions upon which the tree is constructed are the same as used for the FSE 
penetration analysis* That is, the ASW-detects-submarine probability is 
assumed not degraded by prior ASW ship losses and that if a submarine attacks 
the ASW, she is assumed to have sufficient torpedoes to do so* Because of 
the method used in deriving the ASW penetration probability for this analysis, 
parasitic torpedo usage could not be determined. It was observed that two 
mutually exclusive events prevailed, either the submarine penetrated the ASW 
or she was sunk* The method used for computing the probability of these 
events occurring was to first determine the probability the submarine was 
sunk. The desirfed probability of penetration is found by complementation* 

The various paths to "submarine death" during the ASW penetration are 
stylized in Figure To The circled integers, (ijt and the circled arabic A, 



A^, in Figure 7 represent branch points corresponding to those of the tree. 



r to branch point s, indexed by the integer i* From Figures 6 and 7 the 
probabilities for an (I,K) submarine are read off: 



Figure 6, 




represents the event "submarine is 




= PDETU(I,K) X VA(K) 



Paq = PBB(I,K) X PAT(I,K) 
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FIGURE 6 

ARBORESCENCE ILLUSTRATING THE SUBMARINE-ASW SCREEN INTERACTION 
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FIGURE 6 (continued) 
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FIGURE 6 (continued) 
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FIGURE 6 (continued) 
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FIGURE 7 

STYLIZED PATHS TO SUBMARINE DEATH DURING THE SUBMARINE-ASW SCREEN INTERACTION 
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PKBB(I,K) 



pto = (1 - PKBB(I,K)) X PCC(I,i) X (1 - PABB(I,K) + PABB(I,K) x 

□ u 

(PKBBG(I) + (1 - PKBBG(D) x (QTC(177) ) ) ) , i=l, 3 

^01 ^ ^ ■ ^08 " P0X’ ^ 

Pg^ = 1 - PKBB(I.K) 

= (1 - PKINS(D) X PSCA(I,i) X (1 - PAA(I) x (PKSC(I) 

44 

+ (1 - PKSC(D) X (1 - QTC(177)))), i=l, 3 

p^^ = PKINSd) 

P = (1 - PBB(I,K)) X PSUBAT(I) x (1 - PKOSC(I)) x (1 - QTU(309)) 

MX 

f 

p = PBB(I,K) X (1 - PAT(I,K)) X PDETU(I,K) x VA(K) 

+ (1 - PBB(I,K)) X (PSUBATd) x (PK05C(I) x PDETU(I,K) x VM(K) 
+ (1 - PK05C(D) x QTU(389)) + (1 - PSUBAT(I)) x PDETU(I,K))) 



The probability that the submarine penetrates the inner screen, P3(I,K), can 
then be found as 



P3(1,K) = 



1 

1 



Pr (Submarine is sunk attempting A5W penetration) 

Pr (a transit from the first to by all 

possible paths) x Pr (a transit from to first 
by all possible paths) + Pr (a transit from 
to by all possible paths not passing 
through any ) 
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Each of these probabilities can be read from Figure 0 to be the following: 
Pr (a transit from the first to by all possible paths) 

“ P4X * P44 ^^ 4 X P44 ^^ 4 X ^ P44 ^ 4 X^^ 

Pr (a transit from to first by all possible paths) 



Jl 1 2 2 3 3 4 

= Pai Pi4 ^ Pa8^P01 Pi 4 Paa^Psi Pi 4 ^ Paa^Pai Pi 4 * Paa Pai Pi4^^^ 



and 



Pr (a transit from to by all possible paths not passing through 

any ) 



Pa 8 ^Pax Paa^Pax * Paa^Pax Paa Pax^^^ 



Convoy attack phase > It is particularly in the convoy attack that the 
flexibility of the SEALIFT game complicates the expected value analysis. 
Possible realizations of this phase depend heavily upon the chosen war 
option as listed in SEALIFT I j^lj and upon the various torpedo salvo sizes 
chosen for each submarine type. In view of these considerations and the 
problem which prompted this expected value analysis, a particular character- 
ization of the convoy attack phase has been chosen. It is assumed that 

lo one type of torpedo will be used against both cargo and 
screen ships, 

2o the Salvo size against cargo ships is CSALVOd) and that 
CSALVOd) = -J- the initial torpedo supply, 

3. the salvo size against screen ships is A5ALV0(I) and that 
ASALVOd) i 1/6 X CSALVOd ), 

4. in firing one C5ALV0 against cargo ships, the combination 
of CSALVO(I) and the probability of kill per torpedo 
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essentially eliminates the probability of sinking no cargo 
ships, and 

5* having fired two salvos at cargo ships, the submariners 
torpedo supply is exhausted® 

Based on these assumptions. Figure 8 is the simplified tree of that 
part of Event 19 (Submarine versus convoy) of SEALIFT from which the convoy 
attack associated probabilities are derived from the first Markov chain in 
this analysis* For an (I,K) submarine sonar conditions the conditional 
probabilities of interest are 

P41(I,K) = Pr (a submarine shoots one salvo at cargo ships in the 

first Markov chain) 



P42(I,K) 



P5(I,K) 



PD5AU(I,K) 

Pr (a submarine shoots two salvos at cargo ships in the 
first Markov chain) 

1 - PD5AU(I,K) 

Pr (a submarine escapes/ she shoots one salvo at cargo 
ships in the first Markov chain) 

1 - P55AU(I,K) 



and 



P6(I,K) 



Pr (a submarine escapes/ she shoots two salvos at cargo 
ships in the first Markov chain) 



= 1 - PD5AU(I,K) X PS5AU(I,K) 

If the expected number of torpedoes fired by an (1,K) submarine at 
cargo ships in the first Markov chain given the submarine fires one and two 
salvos are E4(I,K) and E5(I,K) respectively, then E4 and E5 are roughly 
approximated by: 



E4(I,K) 



CSALVO(I) 
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FIGURE B 

ARBORESCENCE ILLUSTRATING THE CARGO SHIP ATTACK PHASE 
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and 



E5(I,K) = T0RPS(I,A5W) - CSALVO(I) - E1(I,K) 

The tree associated with the cargo attack phase for those submarines 
which successfully enter the second Markov chain is obtained from Figure 9 
by changing the entire right branch at the branch point K into the event 
“Submarine Escapes*” The conditional probabilities of interest are then 

P4(I,K) = Pr (a submarine shoots her remaining torpedoes at cargo 

ships ) 

= 1 

P11(I,K) = Pr (a submarine escapes given she shoots) 

= P5(I,K) 

The expected number of torpedoes a submarine fires at cargo ships in 
the second Markov chain given she fires, E6(I,K), is roughly approximated by 
E6(I,K) i T0RPS(I,A5W) - CSALVO(I) - 2 x E1(I,K) 

The transition probabilities * The calculation of the Markov chain 
transition probabilities needed for this expected value analysis can now be 
obtained by combining the conditional probabilities derived in the previous 
sections* These transition probabilities are: 

Pi = qi = P1(I,K) 

P 2 = P 2 = (1 - P1(I,K)) X P9(I,K,n) 

= (1 - P1(I,K)) X P7(I,K,n) X P3(I,K) x P41(I,K) x P5(I,K) 

p^ = (1 - P1(I,K)) X P7(I,K»n) X P3(I,K) x (1 - P5(I,K)) 

p^ = (1 - P1(I,K)) X P7(I,K,n) X P3(I,K) x P42(I,K) x P6(I,K) 

Pg = (1 - P1(I,K)) X P7(I,K,n) X P3(I,K) x P42(I,K) x (1 - P6(I,K)) 
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(1 - P1(I,K)) X P7(I,K,n) X (1 - P3(I,K)) 



Pq = = (1 - P1(I,K)) X P0(I,K,n) 

= (1 - P1(I,K)) X P7(I,K,n) x P3(I,K) x P4(I,K) x P5(I,K) 
q^ = (1 - P1(I,K)) X P7(I,K,n) x P3(I,K) x P4(I,K) x (1 - P5(I,K)) 



3p2o3 Expected Number of Torpedoes Fired at Cargo Ships 

In order to determine the expected number of torpedoes fired at cargo 
ships for any given cycle, define the following to be 

F(u,v,I,K) = event an (I,K) submarine fires u salvos at cargo ships 
in the v th Markov chain, u, v =1, 2 
N5(w,I,K) = the number of (I,K) submarines which arrive on station 

for the w th cycle, w=l, 2, 3, poo 
NT(w,I,K) = the total number of torpedoes fired at cargo ships during 
the w th cycle by (I,K) submarines, w = 2, 3, ••• 

The following equation is then used 
(3p20) 

E(NT(w,I,K)) = E(NS(w,I,K)) x |Pr(F(l,l,I ,K) ) x E4(I,K) 

+ Pr(F(2,l,I,K)) X E5(I,K) 

+ Pr(F(l,2,I,K)) X E6(I,K)| 

For the first cycle, distinction is made for the initially on-station, 
initially not-on-station contingencyp For each submarine type and for each 
sonar condition the expected number of torpedoes fired in each contingency 
is computed as 

E (NTj^ (1, 1 ,K ) ) = E (number of torpedoes fired by (I,K) submarine initially 

on station during the first cycle) 
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= E2(I,K) X j PA2(I,K) X E4(I,K) 

+ PA3(I,K) X E5(I,K) 

+ PA5(I,K) X E6(I,K)| 

EtNT^lli I ,K) ) = E (number of torpedoes fired by (I,K) submarines 

initially not-on-station during the first cycle) 

= E3(I,K) X P0(I,K) X I PA7(I,K) x E4(I,K) 

+ PA8(I,K) X E5(I,K) 

+ PAlOd.K) X E6(I,K)j 

The expected number of survivors from the first cycle, ioS* , those who safely 
reach home port, EC(1,I,K), is 

EC(1,I,K) = (E2(I,K) X PA4(I,K) + E3(I,K) x PQ(I,K) x PA9(I,K))x PQ(I,K) 

According to the rules of SEALIFT, the survivors are re-apportioned among 
the sonar boxes to maintain balance® The expected number of re-apportioned 
survivors of the (I,K) classification, EC(2,I,K), is 

EC(2,I,K) = EC(1,I,K) X ^P'I»B0X(J) 

where the sum is over those sonar boxes of like sonar conditions® To deter- 
mine E (NT(2,I ,K) ) , equation (3®20) is applied again, and becomes 

E(NT(2,I,K)) = X EtNT^dd.K)) 

The expected number of survivors from the second cycle, EC(3,I,K), i* 

EC(3,I,K) = EC(2,I,K) x PA9(I,K) x P0(I,K) x P0(I,K) 

For the w th cycle, w =3, the expected number of re-apportioned survivors 
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from the previous cycle, EC(2w-2,I ,K) , serve as the submarine input, and 
the following iteration formulae are found. 



EC(2w-2,I,K) 



EC(2w-3,I,K) 



EC(2w-4.I,K) 

EC(2w-5,I,K) 



and 



E(NT(w,I,K)) 



EC(2w-2,I,K) 



E(NT(w-l,I,K)) 
^ EC(2w-4,I,K) 



The expected number of survivors from the w th cycle (not yet re-apportioned 
among the sonar boxes), EC(2w-l,I,K) is 



EC(2w-l,I,K) 



£C(2w-2,I,K) 



EC(2w-3.I .K) 
EC(2w-4,I,K) 



An estimate of the number of cycles for the war is given by the minimum 
number of cycles a live submarine, not initially on station, might expect to 
achieve® This number is a function of submarine type and is given by the 
integer L(I) 



L(I) 



TWAR 

ALTIM(I) + DTSPRT(I) + 



D5UB5TA 

SUBSPD(I) 



An estimate for total expected number of torpedoes fired at cargo 
ships by (I,K) submarines for the war, E(NT(I,K)), is then 

E(NT(I,K)) = E(NT (1,I,K)) + E(NT_(1,I,K)) + 7 E(NT(w,I,K) 

w=2 

Multiplying this number by the single shot kill probability, PSCS(I,K), pro- 
vides an estimate for the number of cargo ships sunk by (I,K) submarines, 
E(CS(I,K)), 

E(CS(I,K)) = E(NT(I,K)) x PSCS(I,K) 
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Finally, summing E(CS{I,K)) over the indicies I and K provides an estimate 



for the total expected number of cargo ships sunk during the entire war, E(CS) 

3 2 

E(CS) = X Z! E(CS(I,K)). 

I«1 K=1 
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4« Conclusions and Recommendations 



As a model purporting to represent a probabilistic real world situation, 
the SEALIFT game, the analysis presented in this thesis remains largely un- 
tested* The opportunity and data base for testing the model’s reliability, 
accuracy and sensitivity to the proposed approximations are, however, readily 
available* If the SEALIFT game is contemplated as an analytical tool in the 
CYCLOPS III study, the accuracy of this expected value model should be ex- 
plored* An acceptably reliable programmed version of this model should 
effect appreciable savings in computer time* 

For her continued patience, encouragement and good humor, I acknowledge 
the support of my wife in the preparation of this thesis* 
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APPENDIX I 



DETAILED DESCRIPTION OF THE SEALIFT GAME 
The following description is quoted from SEALIFT I j^l 

IIo GENERAL DESCRIPTION 



Convoys 

In this game, friendly forces have two ports to transit between home 
port and delivery port© During this transit the convoy may cross up to two 
“constant” barriers (oneway transit)® The constant barrier in this simulation 
can sink ships that cross it but cannot itself be weakened by these sinkings 
or by the actions of the ships crossing it® Mine fields, air barriers, etc®, 
can be described as “constant” barriers® The “variable” barrier, used against 
submarines, can sink ships that traverse it but can be weakened by having the 
ships that comprise it sunk® Barriers of this type might be submarine bar- 
riers, surface ship barriers, etc® The “variable” type of barrier will be 
described in greater detail later in this segtion® It should be pointed out 
that while each barrier that could be encountered by a convoy or submarine 
will be described herein, any one or all of these barriers may be removed for 
a desired play simply by setting to zero its probability of detecting a pro- 
spective “transitor” • • ® 

At certain times, which are inputs, a convoy is formed in the home port* 
The convoy is made up of three types of ships: cargo ships, ASW escort ships, 

and forward screen ships*® Each convoy when it forms draws its required 
complement of ships from a “port pool®” This pool is stocked by ships from 
returning convoys and by the ship building rate® If the required ships are 
available, the convoy sails toward its delivery port: if not, the game is 

stopped, and a “printout” is made of the game results to date including the 
event that caused the stoppage® 

At one hundred miles from its home port, all ships of the convoy except 
the forward screen escorts pass through a constant barrier (barrier 1C)® The 
losses, if any, are recorded and subtracted from the convoy® 

After having transited the first constant barrier, the convoy becomes 
liable to an air attack® This attack is played as the probability of a con- 
voy ship being sunk per day at sea by air attack® The convoy then enters 
the submarine patrol area and may encounter a submarine while in this area® 

The submarine/convoy encounter will be described in greater detail later® 

Ships lost because of any one of the aforementioned reasons are noted 
by type of ship sunk, method of sinking, and time of loss, and are subtracted 
from the convoy to be printed out at the close of the simulation. If these 



* Variations in inputs permit the forward screen to be composed alternately 
of aircraft, sonobuoys , submarines, etc. 
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losses include ASW escorts or forward screen escorts, the convoy ability 
to detect** enemy submarines will be diminished. Each time ASW escorts 
and forward screen escorts are sunk, new lower detection probabilities for 
the various parts of the convoy are calculated. 

After passing through the submarine patrol area, the convoy becomes 
liable to another air attack, then passes through a second constant barrier 
one hundred miles from its delivery port. After reaching the delivery port, 
the convoy unloads and stays in port a prescribed number of days. During 
this period the cargo ships are again liable to air attack. Of those ships 
sunk in an air attack, one-half are assumed to be sunk before unloading and 
one-half after unloading. The number of cargo ships delivered are now 
tabulated and the convoy is ready to proceed home. 

Proceeding home, the convoy follows the same route in reverse order, 
passing through the same barriers, the same submarine patrol area, and the 
same series of air attacks. When it reaches its home port, the number of 
ships remaining in the convoy are placed in a ”ship pool” by ship type. 

These ship pools will then be used to make up future convoys. 

Submarines 



This model is designed to handle three types of submarines. All of the 
parameters which describe the submarine types are inputs to the game. At 
the start of the war, the opposition submarines, chosen by type of submarine 
as well as by percentage of each type, are placed in two predetermined ”game 
positions,” selected according to the option of the user. The two game po- 
sitions are "on station" and "in port." For those submarines that start the 
game in port, a given number (input) will leave at a given time interval 
(input) until all have left for their stations. 

While a submarine is in port it may be susceptible to strikes at source. 
In addition, the submarine building rate creates new submarines as the war 
progresses. As each new submarine or "turnaround" submarine leaves port, 
it is given a full load of anti-A5W and anti-shipping torpedoes, and a time 
at which it must depart from station and return home. Both of these inputs 
are a function of submarine type. Submarines that are on station at the 
start of the war are given a full load of both types of torpedoes, but their 
time remaining on station is determined by a uniform random distribution. 

The enemy submarine patrol areas are considered to be a constant distance 
from the submarine port. 

While all submarines are considered to transit the same average distance 
to their patrol area, they can appear in one of four areas along the convoy 
route, each area involving either good or bad sonar conditions. Each area 
is looked upon as a "sonar box” and is given a certain percentage of each 
type of submarine (input). A convoy must pass through the sonar boxes in a 



** The terms "detect” or "detection” as used throughout this research contri- 
bution refer to the process of sensing a signal and classifying it as due 
to a target. 
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fixed sequence, viz,, Box 1, first. Box 2, second, etc. Either good or 
bad sonar conditions may be assigned to each box. If a given submarine is 
put into a certain box, the submarine derives both its detection and sinking 
probabilities with respect to the convoy, not only from its type, but also 
from the sonar conditions assigned to the box. Thus for each decision point 
in the submarine/convoy engagement, there is in effect a 2 x 3 probability 
matrix, the rows of which correspond to the two possible sonar conditions, 
and the columns to the three possible types of submarines. 

When the submarine leaves port, it may be required to pass through as 
many as eight barriers, four constant and four variable, before arriving in 
the submarine patrol area. Each barrier can be used or not used as desired 
(input) o The constant barriers are of the same type as those described for 
the convoy, except that in addition to having a probability of being sunk 
as it transits the barrier, the subrftarine can also incur a time delay as it 
goes through or around the barrier. Any time delay so incurred will of 
course cut down the remaining time the submarine can spend in its patrol 
area. 



When the submarine crosses a variable barrier, it not only haS the chance 
of being sunk itself, but also has the chance of sinking a ship in the barrier. 
If the submarine does sink a ship in the barrier, the probability of detection 
for this barrier is assigned a new lower value. The barrier ship sunk is 
replaced at some later time in the game (an input), and the detection prob- 
ability of the barrier increased. In this way the variable barrier not only 
destroys some of the submarines that must transit it, but is itself sub- 
jected to possible reductions in effectiveness. Of course, submarines which 
are detected by the barrier but not sunk suffer time delays which reduce 
the time that can be spent in the patrol area. 

As long as the submarine is at sea, it becomes susceptible to being found 
and engaged by a HUK group. There are two different probabilities per day at 
sea that a submarine will encounter a HUK group; one applies when the sub- 
marine is in transit to or from its patrol area, and the other when the sub- 
marine is in its patrol area waiting for a convoy. Should the submarine be 
found by a HUK group, the HUK group will attempt to sink the submarine. In 
turn, the submarine has a chance of sinking one of the HUK ships. 

The time of arrival ”on station'* depends upon the speed of the submarine 
(which in turn depends upon submarine type) and upon the time delays the 
submarine undergoes while crossing barriers and fighting HUK groups. The 
submarine remains on station or in its patrol area (sonar box) until a con- 
voy crosses the ”subline," The length of the "subline" or what could better 
be called the submarine patrol area, is an input. 

It is assumed that a convoy crosses the center of the submarine patrol 
line. When this happens, the first submarine is drawn from Box 1 and 
placed, using a uniform random distribution, on the subline. Its position 
on this line, relative to the center of the line, is compared to an athwart- 
ship detection range for this type of submarine and for its sonar box. If 
it is within the detection range, the submarine detects the convoy, other- 
wise it does not. If it detects the convoy, it will attempt to "close” on 
this convoy; if it does not detect the convoy, it will be returned to its 
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sonar box (patrol area)* 

The first three boxes are played in this manner: the fourth box is 

played slightly differently* As the convoy approaches the fourth box 
(fourth patrol area), there is given a probability per month of the war 
that submarines in this patrol area will contact the convoy. If contact 
is made, this box will be played as the other three* If, however, this 
contact is not made, this box will not be played for this convoy* This 
fourth box might then be used for a submarine wolf pack; for, although a 
wolf pack has a low probability of intercepting a convoy, given this in- 
terception it has a high probability of sinking ships in the convoy* 

In general, if a submarine has detected a convoy, it will next try 
to penetrate the convoy defenses and sink cargo ships. (However, an option 
is permitted, if desired, by which submarines will try to sink A5W escorts 
in an anti-ASW escort campaign.) The event that plays the submarine en- 
gaging the convoy is Event Subroutine No. 19, "Submarine vs* Convoy." In 
this event, the submarine must first penetrate the forward screen (or out- 
ward screen of the convoy). Given a successful penetration of this screen, 
it must then penetrate the A5W screen (or innermost convoy screen). When 
it has successfully penetrated both of these screens, it may then attack 
cargo ships* When the submarine has sunk one or more ships, either cargo 
or screen ships, the surviving screen ships form surface attack units 
(5AU*s) which attempt to locate the submarine and attack it* Anti-shipping 
and anti-ASW torpedoes are used during this attack, and whenever the allow- 
able minimum number of torpedoes of the type needed to continue the attack 
has been reached, the submarine is required to break off the attack and 
start its transit for home* ..* 

After a submarine has either reached ijbs minimum number of torpedoes 
or has used its allowable time, whichever occurs first, it must return to 
its home port* On the way home, the submarine must traverse, in reverse 
order, the same barriers that it crossed on its way out to station* When 
the submarine arrives at its home port, it will stay there for a certain 
number of days (input) and then be resupplied and again sail for its patrol 
area* 
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